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Abstract. We consider systems of spatial random permutations, where permu- 
tations are weighed according to the point locations. Infinite cycles are present at 
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This article is devoted to random permutations on countable sets that possess a spatial 
structure. Let x be a finite set of points xi,...,xn G M d , and let Sn be the set of 
permutations of N elements. We are interested in probability measures on Sn where 
permutations with long jumps are discouraged. The main example deals with "Gaussian" 
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with [3 a parameter. But we are also interested in more general weights on permutations. 
In addition, we want to allow the distribution of points to be random and to depend on 
the permutations. 

We are mostly interested in the existence and properties of such measures in the ther- 
modynamic limit. That is, assuming that the N points x belong to a cubic box A C 
we consider the limit |A|, N — ► oo, keeping the density p = N/\A\ fixed. The main question 
is the possible occurrence of infinite cycles. As will be seen, infinite cycles occur when the 
density is larger than a critical value. 

Mathematicians and physicists have devoted many efforts to investigating properties 
of non-spatial random permutations, when all permutations carry equal weight. In par- 
ticular, a special emphasis has been put on the study of longest increasing subsequences 
[U |3] and their implications for such diverse areas as random matrices j3j |21] , Gromov- 
Witten theory [22] or polynuclear growth [8], and spectacular results have been obtained. 
The situation is very different for random permutations involving spatial structure; we 
are only aware of the works |16| \T7\ [TO] (and [H]). This lack of attention is odd since 
spatial random permutations are natural and appealing notions in probability theory; it 
becomes even more astonishing when considering that they play an important role in the 
study of quantum bosonic systems: Feynman [9] and Penrose and Onsager [23J pointed 
out the importance of long cycles for Bose-Einstein condensation, and later Siito clarified 
the notion of infinite cycles, also showing that infinite, macroscopic cycles are present in 
the ideal Bose gas |25} 156] . These works, however, never leave the context of quantum 
mechanics. We believe that the time is ripe for introducing a general mathematical frame- 
work of spatial random permutations. The goal of this article is to clarify the setting and 
the open questions, and also to present some results. 

In Section [2j we introduce a model for spatial random permutations in a bounded 
domain A C M. d . As stated above, the intuition is to suppress permutations with large 
jumps. We achieve this by assigning a "one-body energy" of the form — x n u\) to a 

given permutation ir on a finite set x. The one-body potential £ is nonnegative, spherically 
symmetric, and typically monotonically increasing, although we will allow more general 
cases. In addition, we will introduce "many-body potentials" depending on several jumps, 
as well as a weight for the points x. 

As usual, the most interesting mathematical structures will emerge in the thermody- 
namic limit |A|, N — > oo, where N is the number of points in x. The ambitious approach 
is to consider and study the infinite volume limit of probability measures. The limiting 
measure should be a well-defined joint probability measure on countably infinite (but lo- 
cally finite) sets x C and on permutations of x. To establish such a limit seems fairly 
difficult; as an alternative one can settle for constructing an infinite volume measure for 
permutations only, with a fixed x chosen according to some point process. We provide 
a framework for doing so in Section [3j and give a natural criterion for the existence of 
the infinite volume limit (it is a generalisation of the one given in [11] ) . This criterion is 
trivially fulfilled if the interaction prohibits jumps greater than a certain finite distance; 
however, its verification for the physically most interesting cases remains an open problem. 

Another option for taking the thermodynamic limit is to focus on the existence of the 
limiting distribution of one special random variable as | A| , iV — > oo. Motivated by its 
relevance to Bose-Einstein condensation, our choice of random variable is the probability 
of the existence of long cycles; more precisely, we will study the fraction of indices that 
lie in a cycle of macroscopic length. The general intuition is that in situations where 
points are sparse (low density), or where moderately long jumps are strongly discouraged 
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(high temperature), the typical permutation is a small perturbation of the identity map, 
and there are no infinite cycles. In Section [4] we give a criterion (that corresponds to low 
densities, resp. high temperatures) for the absence of infinite cycles. 

On the other hand, infinite cycles are usually present for high density. The density 
where existence of infinite cycles first occurs is called the critical density. We establish the 
occurrence of infinite, macroscopic cycles in Section[5]for the case where only the one-body 
potential is present, and where we average over the point configurations x in a suitable 
way. An especially pleasing aspect of the result is the existence of a simple, exact formula 
for the critical density. It turns out to be nothing else than the critical density of the ideal 
Bose gas, first computed by Einstein in 1925! The experienced physicist may shrug this 
fact off in hindsight. However, it is a priori not apparent why quantum mechanics should 
be useful in understanding this problem, and it is fortunate that much progress has been 
achieved on bosonic systems over the years. Of direct relevance here is Siito's study of the 
ideal gas [26 , and the work of Buffet and Pule on distributions of occupation numbers [6]. 

Section [6] is devoted to the relation between models of spatial random permutations 
and the Feynman-Kac representation of the Bose gas. We are particularly interested in 
the effect of interactions on the Bose-Einstein condensation. While this question has been 
largely left to numericians, experts in path-integral Monte-Carlo methods, we expect that 
weakly interacting bosons can be exactly described by a model of spatial permutations with 
two-body interactions. An interesting open problem is to establish this fact rigorously. 
Numerical simulations of the model of spatial permutations should be rather easy to 
perform, and they should help us to understand the phase transition to a Bose condensate. 

In Section [7] we simplify the interacting model of Section [6} It turns out that the largest 
terms contributing to the interactions between permutation jumps are due to cycles of 
length 2. Retaining this contribution only, we obtain a toy model of interacting random 
permutations that is simple enough to handle, but which allows to explore some of the 
effects of interactions on Bose-Einstein condensation. In particular, we are able to compute 
the critical temperature exactly. It turns out to be higher than the non-interacting one 
and to deviate linearly in the scattering length of the interaction potential. This is in 
qualitative agreement with the findings of the physical community [2j ITS! [15] [20] . In 
addition, we show that infinite cycles occur in our toy model whenever the density is 
sufficiently high. However, our condition on the density is not optimal — higher than the 
critical density of our interacting model. Our condition is not optimal. We expect the 
existence of infinite cycles right down to the interacting critical density, but this is yet 
another open problem. 
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2. The model in finite volume 

Let A be a bounded open domain in and let V denote its volume (Lebesgue measure). 
The state space of our model is the cartesian product 

n AiN = A N x S N , (2.1) 

where Sjy is the symmetric group of permutations of {1,. . . ,N}. The state space Qa,n can 
be equipped with the product cr-algebra of the Borel cr-algebra for A^, and the discrete 
cr-algebra for 5 at. An element (x\, ...,xn) x ir £ JV is viewed as a spatial random 
permutation in the sense that Xj is mapped to x^^j) for all j. Figure |T] illustrates this. 
The probability measure on &a,n is obtained in the usual way of statistical mechanics: 




Figure 1 . Illustration for a random set of points x, and for a permutation 
7T on i. Isolated points are sent onto themselves. 



a reference measure, in our case the product of Lebesgue measure on A^ and uniform 
measure on permutations, is perturbed by a density given by the exponential of a Hamil- 
tonian, i.e. a function H : £l\ t N —> ( — oo, oo]. We will shortly specify the shape of relevant 
Hamiltonians. 

We are interested in properties of permutations rather than positions, and we only 
consider random variables on <Sjv. We consider two different expectations: E x , when 
positions x G A^ are fixed; and -Ea,a> when we average over positions. For this purpose 
we introduce the partition functions 



Y(x) 



Z(A,N) 



7r£iSjv 
1 

M 



-H(x,tt) 



(2.2) 



Y(x)dx. 



A N 



In the last line, dx denotes the Lebesgue measure on Mr . The factor 1/Nl implies 
that Z(A, N) ~ e Vq for large V, N, and for "reasonable" Hamiltonians — a desirable 
property in statistical mechanics. Then, for 9 : Sn -*Ia random variable on the set of 
permutations, we define 



E x (6) 



1 



Y(x 



6(ir)e 



-H(x,tt) 



(2.3) 
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1 



E x (6)Y(x)dx. 



Z(A,N)N\ J a n 

We will be mostly interested in the possible occurrence of long cycles. Thus we introduce 
a random variable that measures the density of points in cycles of length between m and 
n: 

Qm,n(*) = v #{* = h ■ ■ ■ ,N : m < ^(tt) < n). (2.5) 
Here, £i(ir) denotes the length of the cycle that contains i; that is, £i(n) is the smallest 
number n ^ 1 such that n^ n '(i) = i. We also have 

1 N 

Qm,n{v) = yJ2 X ^ m M (^W) ( 2 - 6 ) 
i=l 

with Xi denoting the characteristic function for the interval I. 

We denote by 7Z n the space of random variables on Sn that are invariant under trans- 
positions. That is, 6 £7Zn satisfies 

9(a-\a) = 0(tt) (2.7) 

for all a, n G 5at. Random variables in 1Z at have the useful property that they do not 
depend on the way the set x = {x\, . . . ,xn} is labeled. Instead, they only depend on 
the set x itself, and in that sense are the most natural quantities to study. Notice that 
Qm,n G K-n. 

We now discuss the form of relevant Hamiltonians. H is given by the sum 

H(x, vr) = H<V(x, tt) + Hik H x , + G(x), (2.8) 

k ^ 2 

where the terms satisfy the following properties. Let x = (x±, . . . , xjy). 

• The one-body Hamiltonian has the form 

N 

H^(x,7T) = Y,^-^ {l ))- (2-9) 
i=l 

We suppose that £ is a spherically symmetric function M. d — ► [0, oo], that £(0) = 0, 
and that is integrable. 

• The /c-body term : fiA,Ar — * K can be negative; is has the form 

H^(x,tt) = V(( Xj ,x Aj) ) j&A ). (2.10) 

Ac{l,...,JV} 
\A\=k 

• The function G : ~R dN — > R depends on the points only. It has no effect on the 
expectation E x , but it modifies the expectation E^iv- 

We will discuss in Section [6] the links between spatial random permutations and the 
quantum Bose gas. We will see that the physically relevant terms are £(x) = |x| 2 /4/3, 
that the interactions are two-body (ifW = for k > 2), and that G(x) = 0. From 
the mathematical point of view it is interesting to consider a more general setting. In 
particular, we can restrict the jumps by setting £(x) = oo for |x| bigger than some cutoff 
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distance R. The effect of G is to modify the typical sets of points. We can choose it such 
that Y{x) = 1. We refer to this case as "Poisson", since positions are independent of each 
other, and they are uniformly spread. The point process for the Bose gas is not Poisson, 
however. The fluctuations of the number of points in a subdomain were studied in |18j . 
They were shown to satisfy a large deviation principle with a rate function that is different 
than Poisson's. 



3. The model in infinite volume 

As usual, the most interesting structures emerge in the infinite volume limit V — > oo 
or, more precisely, the thermodynamic limit V — > oo, N = pV. The easiest way to take 



this limit is to consider a fixed random variable, e.g. g ab (ir) from (2.5), and study its 
distribution as V — > oo and N = pV. We will indeed do this in Section |Sj an advantage of 
this approach is that we do not have to worry about infinite volume probability measures. 
But these infinite volume measures are very interesting objects to study directly, in the 
same spirit as when constructing infinite volume Gibbs measures. We advocate this point 
of view in the present section, and introduce a framework for spatial permutations in 
unbounded domains. 

3.1. The cr-algebra. The present and the following subsection contain preparatory re- 
sults about permutations on N. We introduce the "cylinder sets" B^j that consist of all 
permutations where i £ N is sent to j £ N: 

B hJ ={7r£S N :ir(i)=j}. (3.1) 

Let £' denote the collection of finite intersections of cylinder sets and their complements. 
One can check that it is closed under finite intersections, and also that the difference of 
two sets is equal to a finite union of disjoint sets. Such a family of sets is called a semiring 
by probabilists. Semirings are useful because they are easy to build from basic sets, and 
because premeasures on semirings can be extended to measures by the Caratheodory- 
Frechet theorem. Let E be the a- algebra generated by the Bij. 

We start by proving a structural lemma that we shall use when extending finite volume 
measures to an infinite volume one. For A\, A[, . . . , A m , A' m C N, let us define 

Bj}~jr = {vr E 5 N : 7r -1 (i) G A { and n(i) G^Ui^ml. (3.2) 

1" m 

One easily checks that any element of the semiring can be represented by a set of the form 



(3.2). Also, the intersection of two such sets satisfies 



r,A\...A m r-\ r,C\...Cm r>AinCi...A m nC m (o q\ 

A' x ...A' m 1 1 a C' x ...C' m ~ a A'^C' x ...A> m ^C' m ■ 

Lemma 3.1. Let A\, A f 1} A2, Af 2 , . . . be finite subsets of N. // 

B A A }- A A T + 
for all finite m, then lim m _+oo B^}'"^™ ^ 0. 

A f- A m 

It is crucial that both Ai and A\ be finite for all i; counter-examples are easily found 
otherwise. For instance, choose A% = N, Ai = {i — 1} for i > 1, and A[ = {i + 1} for all i; 
then each finite intersection is non-empty. The infinite intersection is empty, on the other 
hand, since there is no possibility left for the preimage of 1. Similarly, choosing A\ = N, 
Aj = {i — 1} for i > 1 and Ai = {i + 1} does not leave a possible image for 1. These two 
cases should be kept in mind when reading the proof. 
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A claim similar to Lemma 3.1 and Theorem 3.2 was proposed in The proof there 



contains a little flaw that is corrected here. 

Proof. We write B®, instead of ^/^'\ ! e * c --- We have 

Bj}-jr = u B a }- a r. (3.4) 

A \-- A m Ol,...,O m O'l-Om V ' 

The union is over ai € Ai, . . . , a' m G A' m , with the restriction that a, / Oj and / a'j for 
i j. The union is disjoint, B a }'" a ] n ^ 0, and 

o 1 ...a m 

B a j~ a r +1 c B , 1 - 8 ™. (3.5) 

a l-" a m+l a l"- a m 

Permutations are charaterised by a "limiting set", namely 

M = (3-6) 

where {a«}, {a'j} are given by 

7r _1 (j) = aj and = a\. (3-7) 
Conversely, given {a^}, {a^}, B a J a J'" is either empty, or it contains the permutation ir that 



satisfies (3.7). We now check that the decomposition (3.4) yields at least one non-empty 
limiting set. 

The sets B®}'"®™ that appear in (3.4) can be organised as a tree. The root is Sj$. The 



sets with m = 1 are connected to Sn, i-e. all the sets B^} , a\ G Ai and G A^. The 

sets with m = 2 are connected to those with m = 1. Precisely, the sets B a J a J , 02 G A 2 

and a 2 G A 2 , are connected to B a j . And so on... There are infinitely many vertices, but 

finitely many vertices at finite distance from the root. 

We can select a limiting set as follows. Let l a y- a j h denote the length of the longest 

a 1 ...a m 

path descending from 5 a ; 1 -- a / m . p or each m, there is at least one {oj},{o^} such that 
gay. .am _ qq_ (Otherwise the tree cannot have infinitely many vertices, since the sets 

a 1 ...a rn 

A\, A[, . . . have finite cardinality.) Further, there exists a m+ i G A m+ i and a' m+1 G A' m+1 
such that £ a j— a y i+1 = OQ< We can choose an infinite descending path such that l a }--- a ™ is 

a l--- a m + l a l--- a m 

always infinite. The sets -B°, 1 "'", m are not empty for any m, so that B a J a J"' is a non-empty 



limiting set. □ 

3.2. An extension theorem. We now give a criterion for a set function \i on £' to 
extend to a measure on E. 



Theorem 3.2. Let Yl be the semiring generated by the cylinder sets B{j in (3.1), and let 
H be an additive set function on S' with fj,(S^) < 00. We assume that for all i G N, and 
all e > 0, there exists a finite set A C N (that depends on i and e) such that 

J2 K B i,j) > 1 - e, and ^ K B j,i) > 1 - £• (3-8) 
jeA jeA 

T/ien /i extends uniquely to a measure on (<Spj,E). 
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If fi is symmetric with respect to the inversion of permutation, i.e. n(Bij) = n(Bj^), 
then the two conditions are equivalent to 

oo 

5>(£y) =m(5n) (3-9) 
j'=i 

for any fixed i. Since 5^ = UjBij, the condition amounts to c-additivity for a restricted 



class of disjoint unions. Thus (3.8) is also a necessary condition. 



Proof. It follows from the assumption of the theorem that for any e > 0, there exist finite 
sets C\, C2, ■ ■ ■ such that, for all i, 

n({ir G 5 N : ir^ii) G Q and vr(i) G C*}) > 1 - 2- i_1 e. (3.10) 

Then for any m, we have 

ri B S>^Z) > 1 " 3* (3-H) 

We prove the following property, that is equivalent to u-additivity in S'. For any decreasing 
sequence (G n ) of sets in £' such that fi(G n ) > e for all n, we have lim n G n 7^ 0. As we 
have observed above, G n can be written as 

G n = B^r-^r (3.12) 

nl ' ' " 71771 

for some finite m that depends on n. Without loss of generality, we can suppose that m ^ n 
(choosing some sets to be N if necessary); we can actually take m = n (by restricting to 
a subsequence if necessary). This allows to alleviate a bit the argument. 



We have G n nG n+ i = G n+ \. Using (3.3), we can choose the sets such that A n i D A n+ \, 



and A' ni D A' n+1 i for each i. The sets A n i,A' ni may be infinite. We therefore define the 



finite sets D ni = A ni D Cj and D' ni = A' ni D Cj. Then 



^MG n )- M ([Bg ::: g : ] c ) (3.13) 



The sets D n i,D' ni are finite and decreasing for fixed i. Let 

A = lim D m , ^ = lim D' ni . (3.14) 

n n 

For any k, there exists n large enough such that 

B Dni-D r c j^Di...Dk c Gfc (3 15) 

The set in the left side is not empty since it has positive measure. Thus the set in the 



middle is not empty for any k. By Lemma 3.1, its limit as k —> 00 is not empty, which 



proves that lirrifc Gk / 0- □ 
Although Theorem |3.2| does not contain any reference to space, we see in the following 



section that the spatial structure provides a natural way for (3.8) to be fulfilled in certain 
cases. 
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3.3. Permutation cycles and probability measure. Here we discuss the connection 
between the considerations of the two previous subsections, and the spatial structure. The 
set of permutations where i belongs to a cycle of length n can be expressed as 

B[ n) = . U . .PI B k _ uji (3.16) 

Jlr-Jn 1=1 



where ji, . . . , j n are distinct integers, and where we set jo = j n . The union is countable 
and j 
then 



In) 

and B- belongs to the a- algebra S. The event where % belongs to an infinite cycle is 



U B\ 

n > 1 1 



(n) 



(3.17) 



and it also belongs to X. We can introduce the random variable £{ for the length of the cycle 
that contains i. It can take the value oo. Since ^ 1 ({ n i) = B^ for all n = 1, 2, . . . , oo, 
we see that it is measurable. 

In general the probability distribution of t% depends on i. Thus we average over points 
in a large domain. Let x C M. d be a countable set with no accumulation points (that is, 
if A C is bounded, then x n A is finite). The elements x±,X2, ■ ■ ■ of x can be ordered 
according to their distance to the origin. More exactly, we suppose that for any cube A 
centered at the origin, there exists N such that x% 6 A iff i ^ N. Let V be the volume of 
A. We introduce the density of points in cycles of length between m and n, by 

d$n(*) = £#{*= 1, ... ,iV : m < 4(ir) < n}. (3.18) 



This expression is of course very similar to Eq. (2.5) for the model in finite volume. 

Next we define the relevant measure on 5^. Let Sjy be the set of permutations that are 
trivial for indices larger than N: 

S N = {tt e S N : 7r(i) =iifi> N}. (3.19) 

We define the finite volume probability of a set B in the semiring £' by 

^\B)=' £ e-"<-* ). (3.20) 

Here, x\ = xf)A. The Hamiltonian H(x\,tt) and the normalisation Y{x\) are given by 
the same expression as in Section [2j 

The existence of the thermodynamic limit turns out to be difficult to establish. If x 
is a lattice such as Z d , or if x is the realisation of a translation invariant point process, 
we expect that v^\b) converges as V — > oo. We cannot prove such a strong statement, 
but it follows from Cantor's diagonal argument that there exists a subsequence (V n ) of 
increasing volumes, such that v^ n \B) converges for all B G £'. (Here, A n is the cube of 
volume V n centered at 0.) Thus we have existence of a limiting set function, v x , but we 
cannot garantee its uniqueness. 

If £ involves a cutoff, i.e. if e - ^) is zero for x large enough, (3.8) is fulfiled and v x 



extends to an infinite volume measure thanks to Theorem 3.2 But we cannot prove the 
criterion of the theorem for more general £, not even the Gaussian. It is certainly true, 
though. 

For relevant choices of point processes and of permutation measures, we expect that 
limy^oo f?m,n(7r) exists for a.e. x and a.e. tt. Let g m n denote the limiting random variable. 
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It allows to define the expectation 

It would be interesting to obtain properties such as concentration of the distribution of 
Q mn . The simplest case should be the Poisson point process. 

3.4. Finite vs infinite volume. The finite volume setting of Section [2] can be rephrased 
in the infinite volume setting as follows. Let fi^' N ^ be the point process on M d such that 



jdx if x C A and I as I = N, 



otherwise. 



Next, let z^j be as in (3.20). For Ai D A2 D A3, let us consider the expectation 

Ea iM ,n(Q&) = J d^ N \x) J dv^\ir)etti>. ( 3 - 23 ) 
This can be compared to Eqs (2.4) and ( |2.5[ ). Namely, we have 

EA, N (em,n) = E aa>n (qW). (3.24) 



It would be interesting to prove that the infinite volume limits in (3.23) can be taken 
separately. Precisely, we expect that 

lim£ A ,p|A| (Qm,n) = A . t 1 ™, A 1 ™. J - E7 Ai,A a ,p|Ai| (0m,n)- ( 3 - 25 ) 

A/R d 1 A 3 /M. d A 2 /'R d Ai/Rd 

4. A REGIME WITHOUT INFINITE CYCLES 

At low density the jumps are very much discouraged and the typical permutations 
resemble the identity permutation, up to few small cycles here and there. In this section 
we give a sufficient condition for the absence of infinite cycles. We consider the infinite 
volume framework. The condition has two parts: (a) a bound on the strength of the 
interactions; (b) in essence, that the points of x lie far apart compared to the decay length 
of e - ^ . Here, B^) denotes the set of permutations where the cycle 7 is present. 

Theorem 4.1. Let x C M. d be a countable set, and let i be an integer. We assume the 
following conditions on interactions and on jump factors, respectively. 

(a) There exists ^ s < 1 such that for all cycles 7 = (ji, . . . ,j n ) C N long enough, 
with ji = i, and all permutations n £ 

n 

Yl \ V (( x ^ x n(j))jeA) - V{{x jlX ^ {j) ) jeA ) < sJ2t( x j- x j-i)> 

AcN j=l 
with xq = x n , and where tt is the permutation obtained from tt by removing 7, i.e. 

\j if j = j k for some k. 



(b) With the same s as in (a) 



71 



e e n 

n ^ ly=(ji,...,j n )k=l 
ji=i 



-(l-*)e(x ik -a:j Jfc _ 1 ) 



< OO. 
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Then we have 



lim lim V u ( x A) (B k} )=0. 

k n 

The condition (a) is stated "for all cycles long enough" , i.e. it must hold for all cycles of 
length n ^ no, for some fixed no that depends on i only. This weakening of the condition 
is useful; many points are involved, they cannot be too close, and m the right side 

cannot be too small. 

If u x extends to a measure on S, then 

lim lim V v£\bW)=v*{B ( F\ (4.1) 

k n 

and the theorem states that i' x (B^°°^) = 0. 

An open problem is to provide sufficient conditions on the Hamiltonian such that, in 
the finite volume framework, we have 

lim lim E K)PV {q m v ) = 0. 



Proof of Theorem 4-1 By definition, 
Vx ^ > - Y(x A ) 2^ ^ 

(4.2) 

with B$ = BW n S N . By restricting the sum over permutations, we get a lower bound 
for the partition function: 

Y(x A )> £ exp{- E ^(fe-^jlW)"^^)}. (4.3) 

Observe that H(x A \y, ir) = H(x\, If). Combining with the two conditions of the theorem, 
we obtain 

E 4 A) 0B? } ) ^ E E f[ e- (1 - s)?( ^-^-i } . (4.4) 

k^n k^n~ l =(j 1) ...j k )l=l 

3l=* 

The right side does not depend on A, and it vanishes in the limit n —* oo. □ 

5. The one-body model 

5.1. Occurrence of infinite cycles. The occurrence of infinite cycles can be proved in 
a large class of models with one-body Hamiltonians, with the critical density being exactly 
known! Here the domain A is the cubic box of size L, volume V = L d . We fix the density 
p of points, that is, we take N = pV. 

Recall the definition (2.5 ) for the random variable g m n that gives the density of points 
in cycles of length between m and n. We consider the Hamiltonian 

N 

H(x,Tr) = E^A( X i ~ X n(i)), (5-1) 
i=l 

where the function £a is a slight modification of the function £, defined by the relation 

e -*A(*) = ^ e-^ x ~ Ly) . (5.2) 
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Our conditions on £ ensure that the latter sum is finite, and that £a converges pointwise 
to £ as V —* oo. This technical modification should not be necessary, but it allows to 



simplify the proof of Theorem 5.1 below. In essence, this amounts to choosing periodic 
boundary conditions. 

Let C = J e" f . We suppose that the Fourier transform of e *W is nonnegative, and 
we denote it Ce" 6 ^' : For k G R d , 

Ce -e(k) = f e -2*ikx e -{(x) dx ( 5 3 ) 



The "dispersion relation" e{k) always satisfies 

• e(0) = 0; e(k) > a\k\ 2 for small k (indeed, the Laplacian of e™ 2 ^ would otherwise 
be zero at k = 0; then J \x\ 2 e~^ x ^ dx = 0, which is absurd). 

• e{k) > uniformly in k away from zero; 

• / e~ E W dk = CT 1 < oo. 

Among many examples of functions that satisfy the conditions above, let us mention 
several important cases. 

(1) The Gaussian: 

e -€(x) = e -N 2 /4/3 ; £ (k) =4ir 2 f3\k\ 2 , C = (4vr/3) d / 2 . (5.4) 

(2) In dimension d = 3, the exponential: 

e -€0») = e-M/^ , e(fe) = 21og[l + (2vr/3|fc|) 2 ], C = 8vr//? 3 . (5.5) 

(3) In <i = 3, a sufficient condition for positive Fourier transform is that 



l( e -t(r)V' 



)" (5-6) 



be monotone decreasing; here, £ depends on r = |x| only. See e.g. |13j and refer- 
ences therein. Thus there exist functions e - ^ with compact support and positive 
Fourier transform. 
(4) In d = 1, 

e -«M = (|x| + I)" 3 / 2 . (5.7) 
Its Fourier transform is positive by Lemma B.l It can be checked that e{k) ~ \k\ 1 ^ 2 



for small k, so that its critical density is finite. 
We define the critical density by 

f dk , , 

J Rd e £{ - k > - 1 

The critical density is finite for d 3, but it can be infinite in d = 1,2. The following the- 
orem claims that infinite cycles are present for p > p c only, and that they are macroscopic. 
It extends results of Stito for the ideal Bose gas [26J, which corresponds to the Gaussian 
case £(x) ~ \x\ 2 . The theorem also applies when p c = 00. 
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Theorem 5.1. Let £ satisfy the assumptions above. Then for all < a < b < 1, and all 

s > 0, 



(a) lim E ApV (g 1Va 

V— >oo 



P ifp^Pc, 
Pc ifp^Pc, 



(b) lim EA, p v(Qv a V*>) = °! 



(c) lim E A:pV (g vb v ) = < 

V— >oo 



ifp^p c ; 

s ifO < s < p - p c , 

P-Pc if0^p-p c ^s. 



The rest of this section is devoted to the proof of this theorem. 



5.2. Fourier representation for spatial permutations. The first step is to reformu- 
late the problem in the Fourier space. In the Gaussian case, = \x\ 2 , this is tradition- 
ally done using the Feynman-Kac formula and unitary transformations of Hilbert spaces. 
There is no Feynman-Kac formula for our more general setting, but we can directly use 
the Fourier transform. This actually simplifies the situation. 

Here and in the sequel, we use the definition (5.3) for the Fourier transform. Let A be 
the unit cube, i.e. we fix the length to L = 1. For / £ L 1 (R d ), let f\{x) = Ylty& d f( x + y)', 
notice that f\ 6 L 1 (A). 

Lemma 5.2. Let f G L 1 (R d ). For all n = 1,2, ... , we have 

,. n 

j dxi . . . dx n Y[ f\{xi - Xi-i) = E f( k ) n - 
jAn i=i fcez" 

(By definition, xq = x n .) 

Proof. The n-th convolution of the function f\ with itself satisfies 



{fT){ x ) = J dxi . . . dx n _i/ A (a;i)/A(x2 - x x ) . . . / A (x - z n _i). (5.9) 
The product J| f\ (xi — Xj_i) is translation invariant, so that 

/ Axi . . .dx„TT/A(xi - Xi-x) = / dxi . . .dx n f A (xi)fA(x2 - Xl)fA(-X n -l) 
jA™ JA n 

= /T(o) 

= E /?(*) (5 " 10) 
= E f^ n - 

k£Z d 

The hat symbol in the third line denotes the L 2 (A)-Fourier transform; but in the fourth 
line, it denotes the L 2 (M rf )-Fourier transform. □ 



We will use a corollary of this lemma. 
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Corollary 5.3. Let A be a cube of size L, A* = be the dual space, and let /a(x) = 

S y ez d f( x + Ly)- Then for any permutation it £ Sn we have: 

. N N N 

/ 6x 1 ...dxtfYif A (xi-x ir ^)= n n^**)- 

•^ AJV i=l ki,...,k N eA* i=l i=l 

Proof. It is enough to consider the case L = 1, the general case can be obtained by scaling. 
The multiple integral factorises according to the cycles of the permutation. Using Lemma 



5.2 for each cycle, we get the result. □ 



We are now in position to reformulate the problem in the Fourier space. By Corollary 



5.3, we have the relation 

1 N 

Ea < n{9) = MX nw £ e{7r) £ e-^-w n (5-n) 

V ' ' n£S N k lt ...,k N eA* 1=1 

with Z'(A,N) = C~ N Z(A,N). 



A simpler expression than (5.11 ) is available for random variables invariant under trans- 
position, i.e. 9 £ TZn- To obtain it, we introduce the set A/a of "occupation numbers" on 
A*. An element n £ A/a is a sequence of integers n = (rife), = 0, 1,2, ... , indexed by 
k £ A*. We denote by Ma,n the set of occupation numbers with total number N: 

M a ,n = [n£ Ma : n k = N }- (5-12) 

fceA* 

Let k = (ki, . . . , fc/v) be an A^-tuple of elements of A*. We can assign an element n = n{k) 
in Na,n by defining, for each k £ A* , 

n k = #{i = l,...,N :ki = k}. (5.13) 

In other words, n& is the number of occurrences of the vector k in k. The map k n{k) 
is onto but not one-to-one; the number of fc's that are sent to a given n is equal to 

M/n nfc! - 

fceA* 

We now introduce probabilities for Fourier modes, permutations, and occupation num- 
bers. The sample space is (A*)^ x Sn', it is discrete, and we consider the discrete a- 
algebra. When taking probabilities, we write k x it for {(fc,7r)}; k for {k} x Sn; n for 
{k : n(k) = n} x Sn and it for (A*)^ x {it}. We define the probability -Pa.at by 

P A ,(h,) = (w e " &£W if = far alii, (514) 
I otherwise. 

Here, = (k%, . . . , fcjy)- Summing over permutations, we get 

= Z'{A^N)N\ 6 ~ Efc6A * £(fc)nfc If. "* ! ( 5 - 15 ) 

with (rife) = n{k). Finally, summing over Fourier modes that are compatible with occu- 
pation numbers yields 

Pa,nH = z^v) e_Efc6A * £(fcK • (5 - 16) 
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It follows that the partition function Z'(A, N) can be expressed using occupation numbers 
as 

Z'(A,N)= ^ e-£*eA* s C0 n *. (5.17) 

neAf A ,N 

These definitions allow to express E\^{6) in an illuminating form. 
Lemma 5.4. For all 9 £ TZn, 

Ea,n(°) = E P A,Jv(n) £ d(TT)P A)N (TT\k), 

with k any N -tuple of Fourier modes that is compatible with n, i.e. such that n{k) = n. 
Proof. Using the definitions ( 5.14[ )-( [5.16[ ), the expectation (5.11) of 6 can be written as 

E K>N {B) = E E e ^) P ^{ k x *0 

= Y P A.^( fc ) E ^) P A,N^\ k )- 



(5.18) 



fce(A*) JV 



The latter sum does not depend on the ordering of the ki's — it depends only on occupation 
numbers (notice that P\ pf(7r\a(k)) = P\^{a~^-K(r\k) for all a G Sn). The lemma follows 
from (5.16). □ 



Lemma 5.5. 



E em,n^) P A,N(^\ k ) = \ E * 



keA* 



n — 171 + I if 1 ^ m ^ n ^ 
rik — m + 1 if 1 ^ m ^ rik ^ n 
if nk < a ^ b. 



Proof. It follows from (5.14) that, given k, permutations are uniformly distributed (over 
compatible permutations). That is, 

PA,N^\ k ) : 



V FlfceA* n k ] - if K(i) = h for all i, 
otherwise. 



(5.19) 



Given k = (ki, . . . , k^r), a permutation tt that leaves it invariant can be decomposed into 
a collection (nk), ^k £ 5 nfc , of permutations for each Fourier mode, namely 

Y Qm,nM P ^MM k ) = ( II E A)ffm,n((Tfe)fc6A*)- ( 5 - 20 ) 



In addition, we have 



Qm^d^keA*) = y E NmA^k) 



(5.21) 



feeA* 



with N m ^ n (iT) being the number of indices that belong to cycles of length between m and 
n. We obtain 



E Sm,n(^) P A,N{A k ) = 7/ E E JV ».»^)' 



(5.22) 
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We see that modes have been decoupled; further, we only need to average N m ^ n (ir) over 
uniform permutations. One easily checks that the probability for an index to belong to a 
cycle of length £, with N indices, is equal to 1/N for any I. Then 



^2 N m,n( n ) =n-m+l, 



for integers 1 ^ m ^ n ^ N . The lemma follows. 



(5.23) 
□ 



We have all elements for the proof of Theorem 5.1 



Proof of Theorem \5.1\ We introduce a set A v of occupation numbers that is typical. With 
po = max(0, p — p c ), let 



A„ = In £ A/a. 



N ■ 



no 
V 



-po 



< rj; n k< vV;n k < V 3v for all \k\ > V'A. (5.24) 

o<|jfc[<y-i ' 

Notice that, for any n G A„, 

i(p, p c )V - 2r]V ^ ^ n k < mm.{p,p c )V + rjV. 



mm 



(5.25) 



It is proved in Proposition A. 2 that -Pa,jv(A)) - *• 1 as V, N — » oo, for all n > 0. Together 
with Lemmas |5.4| and |5.5| we obtain, with N = pV, 



lim E A>pV ( S ) = lim V P AjN (n) Q$ n (n) + Q$ n (n) + Q$ n (n) . (5.26) 



neA v 



We partitioned A* into the disjoint union of = {0}, AfW = {k : < \k\ < V~ v }, and 
M( 2 ) = {k : |fc| ^ V^ 7 ?}, and we introduced 

£m,n( n ) = y [( n - m + 1 ) x [n,oo)("-fc) + (nfe-m + l)X [mi „ ) (n fc ) . (5.27) 

First, we note that Qm]n(n) < V for all m,n, so this term does not contribute in (5.26). 
Next, assuming that 3n < a < b, we have 



min(p, p c ) - 2r/ sC Q{y a {ri) = — ^ n fe < min(p, p c ) + 77. 



(5.28) 



|fe| ^ v-i 



J 1 (2) (2) 

We used (5.25 ). In addition, we observe that Q ya vb (n) and Qyl sV (n) are zero for n € .A^. 
Let us turn to Qm.n- We have 



e W(n) + Q^ vb (n)<V 



(0) 



so these terms do not contribute in (5.26). Finally, we have 



,(0) 

■V b ,sV x 







if p < Pc 



f - if p ^ p c and n < sV, 

s — V y 1 if p ^ p c and no ^ sV 7 ! 



(5.29) 



(5.30) 



Inserting these informations in (5.26) yields the theorem. 



□ 
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6. The quantum Bose gas 

Einstein understood in 1925 that non-interacting bosons undergo a phase transition 
where a single particle state becomes macroscopically occupied. Ever since the Bose 
gas has stirred the interest of physicists and mathematical physicists. The first question 
that theoreticians needed to resolve was whether Bose-Einstein condensation also occurs 
in interacting systems, and what is the order parameter. Feynman introduced in 1953 
what is now refered to as the Feynman-Kac representation [9]. It involves "space-time" 
trajectories and permutations, and Feynman emphasised the role played by long cycles. 
The correct order parameter, called "off-diagonal long-range order", was proposed shortly 
afterwards by Penrose and Onsager [23] . But cycles did not disappear, and it was suggested 
in 1987 that winding cycles are related to superfluidity [23]. Siitd made precise the notion 
of infinite cycles and he showed that it is equivalent to Bose-Einstein condensation in 
the ideal gas [25] [26] . An exact formula made this relation more explicit [27] . Recently, 
Boland and Pule obtained puzzling results for the infinite-hopping Bose-Hubbard model 
[5]: Infinite cycles are present iff there is Bose-Einstein condensation, but the density of 
particles in infinite cycles is strictly greater than that in the condensate. 

There are many open questions, of interest to physicists and mathematicians. The 
relationship between Bose-Einstein condensation, superfluidity, and infinite cycles, has yet 
to be clarified. The ideal gas can condense (if d ^ 3) but it is never superfluid. Interacting 
bosons in one or two dimensions do not condense but they may be superfluid. Infinite 
cycles seem related to Bose-Einstein condensation. It was argued in |28j . however, that 
infinite cycles may be present in a solid at low temperature, so they do not automatically 
imply the existence of a condensate (nor a superfluid). On the other hand, it is expected 
that interacting bosons in three dimensions have a transition to a Bose condensate and 
a superfluid, and that this transition takes place at the same critical temperature. We 
also expect the critical temperature for infinite cycles to be the same. It is therefore of 
physical relevance to consider infinite cycles in the Feynman-Kac representation of the 



Bose gas. We discuss in Subsection 6.1 the relation between the Bose gas and models of 



spatial permutations. The latter should help us understand the effects of interactions on 



the critical temperature. This is discussed in Subsection 6.2 

The Bose gas has also been the object of many interesting mathematical studies. We 
only refer to [30] for a discussion of Bogoliubov theory and related work, and to [19] for 
remarkable results on the ground state of the interacting gas. 

6.1. Feynman-Kac representation of the Bose gas. The state space for a system 
of N identical bosons in a cube A C M. d (size L, volume V = L d ) is the Hilbert space 
L^ ym (A N ) of complex square-integrable symmetric functions on A^. The Hamiltonian is 
the Schrodinger operator 

N 

i=l 1 ^ i,j<N 

Here, Aj denotes the d-dimensional Laplacian for the i-th variable, and U{xi — Xj) is 
a multiplication operator that represents the interaction between particles i and j. We 
choose the self-adjoint extension for periodic boundary conditions. 

In order to get cycles, and following [25], we work in the Hilbert space L 2 (A N ) and we 
consider the unitary representation of Sjy- That is, let U n denote the unitary operator 

Unf{xx, ...,X N )= /(x^!),. . . ^(AT)). (6.2) 
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If 9 is a function of permutations, we can define its expectation by 

<*>A,* = W ^ W] E V, e"^ . (6.3) 

The Feynman-Kac formula expresses the operator e~@ H as an integral operator, whose 
kernel is given by 

K(xi, . . .,xn;Vi, ■ ■ - ,Vn) 

expj- }2 / CMa)-w;(a))dsJ. (6.4) 



1 < i<j ^ N ' 



Here, the Wiener measure Wxy describes a Brownian bridge between x and y in time (3, 
with periodic boundary conditions. We refer to [12] for an excellent introduction to the 
subject. We have 

! dWi y (cj) = £ g p {x -y + Lz)^ gf\x - y), (6.5) 



where gp denotes the normalised Gaussian function 



g (x)= - L-- e-N a /2/3, (6.6) 
^ ' (2vr/?) d / 2 V ; 



The sum over z in (6.5) accounts for periodic boundary conditions. For large L the 
functions (A and gp are almost identical. If / is a function R dn -> E, and < h < ■ ■ ■ < 
t n < (3 are ordered "times" , we have 

/ dW^H/^i),...,^)) = V / da?i...dx n 

J ' z6Z-- /R<,n 

% (s - xi)g^l tl (x 2 - xi) . --gfl^y - x n )f(x u . . .,x n ). (6.7) 
By the Feynman-Kac formula, the partition function is 

Z(A, N) = Tr L| ym (A^) e_/3H = XI J^dxi...dx N K(xi,...,x N ;x^,...,x^). 



1T<=S 



(6.8) 

The sum over permutations is present because we work in the symmetric subspace. The 
expectation of observables on permutations can be expressed as 

(0)a,n = z , k ^2 j^dxi...dx N K(xi,...,x N ]Xi T{ i ) ,...,x 7T ( yN) ). (6.9) 

The connection to the model of spatial permutations is immediate in absence of inter- 
actions (U(x) = 0): (6) a ,n = E A)N (9) with £(jc) = \x\ 2 /4(3, ff( fc ) = for k > 2, and 
6?(as) = 0. 
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6.2. Discussion: Relevant interactions for spatial permutations. The models of 
spatial random permutations retain some of the features of the quantum Bose gas in the 
Feynman-Kac representation, but not all of them. The interactions between quantum 
particles translate into many-body interactions for permutations. However, an expansion 
reveals that to lowest order the interaction is two-body; precisely, the interaction between 
jumps x i — ^ y and x' *— > y' is given by 

V(x,y,x',y') = J [I - e-if^-W^Jd^^^w) (6.10) 

for x ^ y. Here, W£ y = 5 < ~ 1 (x — y)W^ y is a normalised Wiener measure, J dW* y = 1. 
If U consists of a hard-core potential of radius a, we notice that V(x, y, x', y') is equal to 
the probability that a Brownian bridge, from x — x' to y — y' , intersects the ball of radius 
a centered at 0. 



The computation of (6.10) can be found in [29]; it is partly justified by a cluster expan- 
sion, although a rigorous result is still lacking. We expect that the critical temperature 
for the occurrence of infinite cycles is the same for {-)a,n and Ea,n(~) to lowest order in 
the strength of the interaction potential. 

An important question about Bose systems concerns the effect of interactions on the 
critical temperature. Over the years physicists gave several, conflicting answers. But a 
consensus has recently emerged in the physics literature, mostly from path-integral Monte- 
Carlo numerical studies. The critical temperature Tc , as a function of the "scattering 
length" a of the interaction potential (see e.g. [19] for the definition), behaves in three 
dimensions as 



rp(a) _ rp(0) 

C T<°> " 



cp^a + oip^a), (6.11) 



with c ~ 1.3. See [TJ] [151 ED] and references therein. The model of spatial random 
permutations should give the exact correction to the critical temperature, i.e. the correct 



constant c in (6.11). 



It is worth discussing the effects of interactions on spatial permutations in detail. They 
both modify the point process and the measure on permutations. Let //W denote the 
point process for an interaction potential of scattering length a, and the measure on 
permutations. We can assign different parameters to these two measures, so as to have 
the expectation 

E ^,a')(0) = [ dp( a \x) [ d4°V)0(7r). (6.12) 



Let us fix the particle density, and let Tc a,a ^ denote the critical temperature for the 
occurrence of infinite cycles. Of course, T c (a) =T c (a ' a) . For small a we should have 

rp(a) _ rp(0) rp(a,0) _ rp(0) rp(Q,a) _ rp(0) 

J-c J-c J c i c -*c c i / \ /„ 1 „n 

go) = — — + — — + ( 6 - 13 ) 

The first question is whether both corrections are linear in o; if not, we could happily 
dismiss one. This formula may be useful in numerical simulations. The second term in 
the right side involves the point process for the ideal gas, while permutation jumps are 
interacting. It can be easily calculated numerically since the set of permutations is discrete. 
To know the corresponding linear behaviour would already provide some information. The 
first term involves the positions of particles, and it may be difficult to generate a typical 



realisation of the points. In any case, the random permutation approach with Eq. (6.10) 
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is much easier than path-integral Monte-Carlo simulations, and it should yield the same 
result to lowest order. 

7. A SIMPLE MODEL OF SPATIAL RANDOM PERMUTATIONS WITH INTERACTIONS 

In the preceding section we have discussed a two-body interaction that is exactly related 
to the quantum Bose gas. Rigorous results seem difficult to get, though. In this section we 
simplify the interaction so as to retain only the largest contribution. Our approximation 
is not exact, but the model may be of interest as an effective model, and it is exactly 
solvable. 

7.1. Approximation and definition of the model. It helps to think of particles as 
describing Brownian motions, and to interact whenever they cross each other's paths. 
Clearly, particles that belong to the same 2-cycle interact a lot. Our approximation 
consists in retaining only these interactions, and in neglecting the rest. Thus we consider 
the Hamiltonian 

1 N 

H(x,ir) = —^2\xi- x^^l 2 + ^ V(xi,Xj,Xj,Xi). (7.1) 

P i=l 1 < i<j < N 



The potential V(-) is the exact jump interaction given in (6.10). To lowest order in the 
scattering length, we find that [1] 

V(x,y,y,x) = r 2 ^ 1 + 0(a 2 ). (7.2) 
\x - y\ 

If 9 is a random variable that depends only on permutations, its expectation is given by 

7T£iSjV 



This allows to simplify the model (7.1 ) further, without additional approximation. Namely, 
we introduce the simpler Hamiltonian 

1 N 

H( a \x, vr) = — \*i ~ **(0l 2 + <*N 2 {ir), (7.4) 



i=l 



with ^(vr) denoting the number of 2-cycles in the permutation ir. Expectations with H 
and are identical provided 

f dxe-^ri = f dxe -H(»Hw,*) (75) 
Ja n Ja n 
for any fixed permutation it. Both sides factorise according to the cycles of tt. The 
contribution of cycles of length 1,3,4,. . . is identical. We then obtain an equation for cycles 
of length 2, namely 

f d Xl dx 2e -^ Xl - X2? - V{x ^ X2 ^ Xl) = [ d Xl dx 2 e-^ lxi - X2l2 - a . (7.6) 
J A 2 Ja 2 



Using (7.2), a few computations give Hj 

1/2 



«= (^J ' a + 0(a 2 ). (7.7) 
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We emphasise that the approximation consists in retaining only interactions within 2- 
cycles. Computations are exact afterwards, at least to lowest order in the scattering length 
of the interaction potential U between the quantum particles. 

7.2. Pressure and critical density. We now generalise a bit the model above by con- 
sidering more general one-body terms, as we have done throughout this article. Let 

N 

H^(x,n) = J2^A(xi-x 7ri{) ) + aN 2 (n), (7.8) 

i=l 

with £a as in Section [5] and a is a positive parameter. 

The pressure of this simple model can be computed exactly. One gets the critical density 
by analogy with the ideal gas. We state below a result about infinite cycles for a larger 



density, see Theorem 7.2. It remains an open problem to show that infinite cycles are 
present all the way to the critical density. 

For given a, the pressure depends on the chemical potential \i and is given by 

P (o) W = Jim ~logZ(A,/i) (7.9) 

V— >oo V 

with Z(A,fi) the "grand-canonical partition function". We can define it directly in the 
Fourier space, by 

uN N 

= E E E e " a " 2W II *- e(h) s ^ y (7-10) 

N ^ fci,...,fcjveA* TreSiv i=1 

We need the pressure of the ideal gas 



P 



(0)/ 



(p) = - / log(l - e~( £ W~^ )dk. (7.11) 



Theorem 7.1. For \i < 0, the limit (7.9) exists, and 



p(«)( M ) = P (°)( M ) - i e 2 ^ [1 - e" a ] / e- 2e(fe) dA: 



Notice that the model is defined only for \i < 0, like the ideal gas. Its derivative at 
/i = 0— gives the critical density, and we find 

>) = ^ 

The first term of the right side is equal to the critical density of the ideal gas, Eq. ( |5.8| ) 
The second term is the correction due to our simple interaction. 



(0) _[!_ e -«] /" e - 2£ ( fc )dfc. (7.12) 



The physically relevant situation is d = 3, £(x) = |x| /4/3, and a in (7.7). We find that, 
to lowest order in a, 



rp(a) _ rp(0) 

C r<°> ' 



gp 1/3 a, (7.13) 



with c = 0.37. The details of the computations can be found in [4J. This formula can be 



compared with (6.11). Our constant has the correct sign, and is about a quarter of the 
expected constant. This suggests that the simplified model accounts for some of the effects 
of interactions on the critical temperature for Bose-Einstein condensation. Besides, it 
allows for simple but illuminating heuristics: Interactions discourage small cycles; then all 
other cycles are favoured, including infinite cycles. As a consequence, repulsive interactions 
increase the critical temperature (or decrease the critical density). 
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Proof of Theorem 7.1 From ( |7.10[ ), we have 



z(a, M )= yi n 

(nfc)iteA* fceA * 



e -(e(fc)-^)n fc J_ e -aiV2(7r fe ) 



(7.14) 



We decomposed the permutation ir into permutations (tt^) for each Fourier mode, and we 
also used 

iV 2 (7r) = N 2^k)- (7.15) 

keA* 

Notice that the chemical potential needs to be strictly negative, as in the ideal gas. We 
get 



p(«)( M ) = lim I y log IV e -( e w-M)n y i 

fceA* L n ^ 7r65 n 



3 -aAf 2 (7r) 



(7.16) 



Let us compute the bracket above. For given 7r S <S n , let denote the number of cycles 
of length j. Then Yljj r j = n i an d the number of permutations for given (rj) is equal to 



The bracket in ( |7.16| ) is then equal to 

El \ - n 

n > ri.ro.... > 1 1 3 * 1 



n ^ ri,r2,... ^ 



E II ^[Je- iW * ) - |t, ] rj E ^[K 2(£(fc) -^T 2 



>"1 i r 3 , r 4 ) ■ • • 5^ 0jr' = l,3,4,. 



r 2 ^0 



= exp{ £ 4 e -mk)-K) + 1 e -2( S (fc)-^)-a | 
3=1,3,4,... 

= exp{-log(l- e -(^)-^)) _ i e -2( £ (fc)-M) [!_ e -«]j. 

We can insert this into (7.16). In the limit V — > oo the expression converges to a Riemann 
integral. □ 

7.3. Occurrence of infinite cycles. Given a £ [0, oo], the expectation of a random 
variable on permutations is given by 



N 



-^A^i-^fi)) _ (7.17) 



The normalisation Z(A, N) depends on a, although the notation does not make it explicit. 
We expect that the claims of Theorem 5.1 extend to a ^ 0, with the critical density given 
by (7.12) instead of (5.8). But we only state and prove a weaker claim. 

Theorem 7.2. For all < b < 1, 



lim E A>pV (Qyb y) ^ p- — 
V— >oo 1 



(1+ e-«) 



2 Pc • 
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Of course, the theorem is useful only if the right side is strictly positive. The proof is 
similar to Theorem 5.1, but there is one important difference. We cannot invoke a set of 
typical occupation numbers, such as A v in Eq. (5.24). We prove below (Proposition 7.6) 
that the zero Fourier mode is macroscopically occupied if the density is large enough. But 
we need it to be strictly bigger than p^ , itself bigger than p^ . 



Let us define 



-aN 2 (ir) 



(7.18) 



TreSn 



Notice that h n (0) = 1. We introduce a probability on (A*) N x Sn that generalises Eq. 

1 



P AtN {k x tt) 



(7.19) 



Z'(A,N)N\ 

if ki = k^u\ for all i, it is zero otherwise. Summing over permutations, and over vectors k 
that are compatibles with n, we get 

Pa > n{ti) = wrm n e ~ e{k) ( 7 - 2 °) 

We generalise Lemma |5.4| 
Lemma 7.3. We have Z(A, N) = C N Z'(A, N) with C = J e~* ; and for all 6 G TZ N , 
EaM ) = Y P A,Jv(n) Y 0(K)P AtN (TT\k) 

with k any N-tuple such that n(k) = n. 



Proof. We use Corollary 5.3 to rewrite the expectation E A jv i n the Fourier space: 

i N 
Ea, n (9) " " 



1 



Z(A,N)N\ 



ttGSn fci,...,fcjvSA* i=l 

= Y 9 W Y P ^ k x vr) 

= Y P *M*) Y G(TT)P A ,N(Mk). 

fce(A*) JV *-eS N 
One can sum first over n and then over compatible fc's. 



(7.21) 



□ 



Next we gather some information on the functions h n (a) defined in (7.18). 
Lemma 7.4. For a G [0, oo], fei 8 = \{1 - e~ a ) 6 [0, §] . 

(a) M«) = Ei!ijr(-' y ) i - 

(b) < /i n (a) < 1. 

(c) e- 5 -<5"/ 2 /L§J! < hn(a) < e~ 5 + 5 n / 2 /L§J!- 

Proof. Isolating the contribution of the cycle that contains 1, we get the following recursive 
relation; for n 3s 2, 



1 n ~ 1 1 
n (a) = - y~] hj(a) (1 - e~ a )/i n _ 2 (a) 



(7.22) 



3=0 
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We also have ho(a) = /11(a) = 1. Now the formula in (a) can be proved by induction, (b) 
is a consequence of the alternating series in (a). Notice that /12(a) = /13(a) = 1 — 5. (c) 
follows from the expression 

h n (a) = e~ 5 - Yl ^Y~- ( 7 - 23 ) 

Recall that alternating series of decreasing terms are bounded by their first term. □ 
Let us define N a ^(ir) = J27=i N °te that iY ai „(7r) = Vq^tt). 



Lemma 7.5. Suppose a > 2, and let 5 as in Lemma \l .J\ For all m ^ 0, 

e -aN 2 (n) 

h n (a)n\ 



E - a N 2 (ir) 
N a , n (7v) > (n - o - 2m)(l - <5 m )- 



Proof. When summing over permutations, all indexes i in the definition of N an (Tr) are 
equivalent, so that 

Summing over the lengths of the cycle that contains 1 , we get 

-aN 2 {n) n p ~aN 2 (n) 



^" w w =E»(- 1 )-(»-i+ 1 ) £ Ma)nl 

J— " 7TfcO n _j 



(7.25) 



We get a lower bound by summing up to n — 2m. By Lemma 7.4 (c), we have for 
2 < j < n — 2m 

hn-j(a) e- s -8™/m\ >1 _ §m (7 2Q) 

□ 



The last result that is needed for the proof of Theorem 7.2 is that our interacting 
model displays Bose-Einstein condensation, in the sense that the zero Fourier mode is 
macroscopically occupied. 

Proposition 7.6. The expectation for the occupation of the zero Fourier mode is bounded 
below by 

Proof. We proceed as in Appendix B of [27]. We have 

E A , N {f)=p-^ Yl E ^M^k), (7.27) 
fceA*\{o} 
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and 



E E < 

j > 1 n£Af A ,N-j 



(7.28) 



Z'(A,N) 



hn k {a) 



The latter ratio is smaller than (1 — 5) 1 by Lemma 7.4 (b). By restricting occupation 
numbers to no ^ j, we also have 



Z f (A,N)> II e ~ £(fcK M 
>Z'(A,N-j)(l-5). 



hn +j(a) 
h no (a) 



(7.29) 



Then 



It follows that 



E A , N (n k ) < (1 - 5)- 2 J] e 



(1-5)- 



1 



3 e(fe) _ 1 ' 



We get the proposition by letting V — > oo. 



keA*\{0} 



Proof of Theorem 7.2. From Lemma 7.3, we have 



Ek,N{Q m ,n) = E P KN{n) ^ Qm^) 

k i= k Ty(i) Vi 



• ttA : 2(7r) 



FUeA* V(«K ! ' 



(7.30) 

(7.31) 
□ 

(7.32) 



Compatible permutations factorise according to Fourier modes, i.e. ir = (tt^) with 7Tfc 6 
5 nfe . Also, iY 2 (vr) = J2 k ^2(vr fe ). Then 



E ^(n)(n E 



-aN 2 (-K k ) 



keA* w k es nk 



h nk (a)n k \ 



Qm,n ((^fc)) 
e -«JVa(irj.) 



(7.33) 



fceA* n k es n . 



K k (a)n k \ 



We keep only the term = 0. Using Lemma 7.5 we obtain the lower bound 



Ea,n{Qv\n)> E ^A,Ar(n) n ° J V (1-S vb ) 



>E AiN (^)-AV b - 1 . 
The claim follows from Proposition |7.6[ 



V 



(7.34) 



□ 
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Appendix A. Macroscopic occupation of the zero Fourier mode. 



In this appendix we investigate the random variable nwno under the measure (5.16) 
in the thermodynamic limit V — > oo, N = pV, for all density parameters p. We want to 
show that uq/V approaches a limit for each density p. We will actually show much more, 
by giving the limiting moment generating function of Hq/V. We partly follow Buffet and 
Pule [6], who considered the ideal Bose gas in arbitrary domains. 



Theorem A.l. Let po = max(0, p — p c ), with p c the critical density defined in (5.8). Then 



lim E\ pV (e 

V^oo 



An /V \ _ Ap 



for all A ^ 0. 



Our first proof applies only when p c is finite. We give below an argument that completes 
the proof. 

Proof when p c < oo. It is shown in [27], Appendix B, that for all k £ A*, we have 

P ^ nk>j) = e -^^EzA. (AJ , 

Since P(rik = j) = P(rik ^ j) — P{nk ^ j + 1), we find for b > 

1 N 

1 ' ' j=o 



^e-^(Z'(A,j)-Z'(A,j-l)). (A.2) 



Z'(A,JV) . „ 

Here, we used the convention Z(A, — 1) := and the fact that -PA,Ar( n fc ^ A + 1) = 0. 
Putting in A = [pVj and setting v = A/V we obtain 

E KpV (e A ^) = z^pV) P Q {Z ' {Kj) ~ Z ' {Kj " (A - 3) 

Above, we wrote pV instead of [pV\ and we will continue to do so, to simplify notation. 
Now comes the clever insight of Buffet and Pule [6 : In (A. 3), both Z'(A, pV) and the sum 
over j can be written as integrals with respect to a purely atomic, ^-dependent measure 
p,\ on IR + ; since the functions that are being integrated will not depend on V, we only 
need to study the limit of p,\. The measure pL\ is given by 

oo 

MA := C A J2(Z'(A,j) - Z'(A,j - l))S j/v 

on M + . S x denotes the Dirac measure at x, and the constant C\ will be fixed later in order 
to obtain a limit measure. From (A. 3) it is now immediate that 

Ek v ( e Xn 0/ V \ _ e Ap fl M (x)e-**^(dx) 



SPATIAL RANDOM PERMUTATIONS AND INFINITE CYCLES 27 

What makes the idea work is that we can actually calculate the Laplace transform of ^a 
and take the limit. We have 

/ e- x ^ A (dx) = C A J2^(Z'(A,j)-Z'(A,j-l)) = 
Jo 



3=0 




00 

j=0 



e v j exp 



* (-Y> g (l-e^ 

V feeA* 



y-e(k) 



fceA*\{o} 

The second equality above is just an index shift, the third is the well-known formula for 
the pressure of the ideal Bose gas, Eq. (7.11 ), and the last line follows from e(0) = 0. The 
sum in the exponent of the last line is actually quite manageable. By the fundamental 
theorem of calculus we have for each k ^ 

1 r x 1 



log (1 - ) = $J „_^ y + ^ (1 - 



and summation over k gives 



Y, ^g (l - e -v-^) ) = Y ^g (l - e^ £ ( fc ) ) + A Pc ,A (A.5) 
fceA*\{o} fceA*\{o} 



with 

Pc,A 



A 



J Q V Yl e y/V+e(k) _ 



fceA*\{o} 



The first term in (A.5) diverges as V — > 00 and defines Ca- The second term converges to 
the critical density p c : the integrand is decreasing as a function of y and converges to p c 
as a Riemann sum for each fixed y, since 



e -y/v 



V ^ e £ ( fc ) - 1 ^ V ^ e y/v+e(k) _ 1 ^ v Z> e e(k) _ 1 ■ 

fceA*\{0} fcGA*\{0} fceA*\{0} 



Dominated convergence in y now proves convergence to p c . We have thus shown that for 
all A > 



poo 

lim / e- Xx fi A {dx) = e- xpc 
Jo 



and thus by the general theory of Laplace transforms p A — > 5 Pc weakly. When used in 
(A.4), this shows the claim for p > p c . For p < p c , both denominator and numerator go 
to zero, and we need a different argument: we note that by what we have just proved 

lim lim E^ pV {e Xno/V ) = 1. 
p\p c V^oo 

Since the expectation above can never be less than one, all we need to show is monotonicity 
in p, i.e. 

P\,N+i(n > j) > PA,N(na > j)- (A.6) 
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For this we use (A.l) and we obtain 



P\,N+i(n k ^ j) = P^ N {n k ^ j) 



P AjN+1 (n k > 1) 
^A,JV-i+l("fc ^1) 



so it will be enough to show (A. 6) for j = 1. By (A.l) this means we have to show 

Z'(A,N) 2 ^ Z'(A,N- 1)Z'(A,N + 1). 

Davies [7] showed that the finite volume free energy is convex, which proves the inequality 
above. □ 

Proof of Theorem \A.l\ when p c = oo. We get from Eq. A. 2 after some rearrangements of 
the terms 



N 



Wv \ _ e -A/v 



(1 _ e - VV ) e V ^ e~ Xi/v e - y[qA{i/v) - qA(N/v)] . (A.7) 

i=0 

q A (p) is convex and its limit q(p) is strictly decreasing for p < p c (and p c = oo here). 
Beside, we have q A {i/V) — q A (N/V) ^ 6 > 0, for all ^ i ^ iV/2, uniformly in V and 
iV = pV . The right side of ( |A.7 ) is then less than 



(1-e 



iV/2 Af 



i=0 



i=N/2 



This clearly vanishes in the limit V — > oo. 



□ 



Recall the definition of the typical set of occupation numbers A„ in (5.24). 
Proposition A. 2. For any density p, and any 7] > 0, 

lim P\ pV (A n ) = 1. 

V— s-oo 

Proof. Let us introduce the following sets of unlikely occupation numbers: 

= (K) : |n /y-po| >r/}, 

^4 (2) = {("*) = n ^ 
o<|jfe[<v-f 

A® = {(n fc ) : n fc > V 3,? for some |Jfc| > y - '}. 



(Ai 



(A.9) 



A.l 



a c v = aWua^ua^. 

that -Pa,pv(^4^) vanishes in the limit V — ► oo. Equation 

:■</-■)' p-V(g A (^)-gA(f )) . ( A . 10 ) 



Then 

It follows from Theorem 
(A.l) can be written using free energies as 

Pa.tvK >i)= e~ £ W l e 

We have already mentioned that q A is convex, and the limit q is decreasing. Then 

qA(p ~ e) - q A (p) ^ u, (A.ll) 

with v 0. In addition, ^ > below the critical density, or if the critical density is 
infinite. Since 

E A ,N(n k ) = PA,N( n k > «)> (A- 12) 

i > 
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we get a bound for the expectation of occupation numbers for k 7^ 0, namely 

E A>N (n k ) < e<k) l_ r (A.13) 



From Markov's inequality and the bound above, we get 
Pa,n(A^) ^ 



Ylo<\k\<v-i E A,N{n k ) 



rjV 



(A.14) 



e e(fc)+y _ x ' 



The right side vanishes as V — > 00: Indeed, this holds if the critical density is finite even 
if v = 0; and we know that v > if the critical density is infinite. 

Finally, (A. 10) implies that P/v,N{ n k ^ *) ^ e — . Since e(A;) > a|/c| 2 for small fc, we 
have for large V 

Pa,jv(^ (3) )< Yl P KN{n k >V^) 

\k\ > V-i 

^ 2^ e (A.15) 

|fe|>V-i 

fceA* 

The prefactor of the last line goes to 0, while the sum converges to a finite Riemann 
integral. The whole expression vanishes in the limit. □ 

Appendix B. Convexity and Fourier positivity 

The case d = 1 of the following lemma was used to provide an example of a one- 
dimensional system with finite critical density. The result is certainly known — we were 
told that it may go back to Polya — but it is easier to find a proof than a reference. 

Lemma B.l. Let g : (0,oo) — > (0, 00) such that f^°r d ~ 1 g(r)dr < 00, and such that 
r d ~ 1 g(r) is convex. Then the function f(x) = g(\x\) on lBL d has positive Fourier transform. 

Proof. First, we show that for any convex function h on (0, 00), we have 



rn+l 

/ h(u) cos(2vru)du ^ (B.l) 

J n 



for any integer n ^ 0. It is enough to consider the case n = 0; the general case follows 
from a change of variables (translates of convex functions are convex) . We have 

1 1-1/4. 

h(u) cos(27ra)du = / [h(u) - h(\ - u) - h(\ + u) + h(l - u)} eos(27m)dw. (B.2) 
Jo 

We used the fact that 

cos(2vra) = - cos(2vr(i -«)) = - cos(2tt(± + u)) = cos(2vr(l - it)). (B.3) 

\ k 1/2 

We now show that the bracket in (B.2) is positive. Let a = 1 _ 2u - Because h is convex, 

/i(tj — u) = h(au + (1 — a)(l — u)) ^ ah(u) + (1 — a)h(l — u), 
h{\ +u) = h((l - a)u + a(l - u)) < (1 - a)h(u) + ah(l - u). 
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This proves positivity of the bracket of ( |B.2 ), hence of (B.2) and (B.l). The case d = 1 



of Lemma 



B.l 



follows, since f(k) = 2 f(r) cos(2-7r|A;|r)dr. 
For d ^ 2, let 9 denote the angle between x and k, and let ^ denote the measure 
of all remaining angles = 2 for d = 2, = 27rsin# for d = 3). Then 



/(*) 



r d-l dr 



t/ 2 - 



l(6)d.Bg(r) cos(2vr|fc|rcos6») 



5(6>)dfl 

pfcos 5) £ 



duu d_:L 



7/ 



(B.4) 



cos 



COs(27Tu). 



Now u d ~ 1 g( 



\k\ cos / 



is convex in u for given \k\ and given < 6 < \ (scaled convex 



functions are convex). The latter integral is positive by (B.l). 



□ 
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